Relaxation in quantum dots due to evanescent-wave Johnson noise 
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We present our study of decoherence in charge (spin) qubits due to evanescent-wave Johnson 
noise (EWJN) in a laterally coupled double quantum dot (single quantum dot). The high density 
of evanescent modes in the vicinity of metallic gates causes energy relaxation and a loss of phase 
coherence of electrons trapped in quantum dots. We derive expressions for the resultant energy 
relaxation rates of charge and spin qubits in a variety of dot geometries, and EWJN is shown to be 
a dominant source of decoherence for spin qubits held at low magnetic fields. Previous studies in this 
field approximated the charge or spin qubit as a point dipole. Ignoring the finite size of the quantum 
dot in this way leads to a spurious divergence in the relaxation rate as the qubit approaches the 
metal. Our approach goes beyond the dipole approximation and remedies this unphysical divergence 
by taking into account the finite size of the quantum dot. Additionally, we derive an enhancement 
of EWJN that occurs outside a thin metallic film, relative to the field surrounding a conducting 
half-space. 

PACS numbers: 03.67.-a, 03.65.Yz, 42.50.Lc, 73.21.-b 
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I. INTRODUCTION 

Semiconducting quantum dots are promising candi- 
dates for scalable quantum information processing^. Sev- 
eral experiments performed on laterally coupled double 
quantum dots (DQDs) have demonstrated precise and 
rapid control of the coupling between electronic charge 
states and coherent manipulation of trapped electrons?^^— 
leading to realization of a DQD as a qubit. Quantum dots 
are realized in a variety of experimental setups, including 
a Si and GaAs two-dimensional electron ga.Sr-'— semicon- 
ductor nanowircs^ and carbon nanotubcsi^ In almost all 
of these implementations, confinement and manipulation 
of an electron in a quantum dot is achieved by applying 
an electrostatic potential through metallic gates. While 
the metallic gates are crucial for qubit control, they can 
also act as a source of decoherence during qubit opera- 
tions. 

Several decoherence mechanisms, such as hyperfine 
coupling of the trapped electron spin to host lattice nu- 
clear spins in spin-based qubitsi and electron coupling to 
phonon modes^"— in charge-based qubits, have been pre- 
viously studied in an effort to identify the major source 
of decoherence in semiconductor qubits. A more recent 
study investigated decoherence due to voltage fluctua- 
tions in the metallic gates using the lumped circuit model 
of a DQD charge qubit In almost all of these stud- 
ies;^""— the estimated energy relaxation rate is at least 
an order of magnitude smaller than the rate observed ex- 
perimentallyf^i^ suggesting that a different decoherence 
mechanism is dominant in current experimental setups 
for charge qubits. 

Here we present our study of decoherence in a quan- 
tum dot due to electromagnetic field fluctuations near the 
metallic gates. We focus primarily on noise from the high 
density of evanescent modes in the vicinity of metallic 
gates. This evanescent- wave Johnson noise (EWJN) has 
been identified as an important source of decoherence in 



atomioi^ii^ and quantum dot based qubitsJ^ Our previ- 
ous worki^ as well as other theoretical estimatesi^ of the 
effect of Johnson noise in atomic and quantum dot based 
qubits use the dipole approximation, which is a valid ap- 
proximation if the distance from the metallic gate to the 
qubit is much larger than the size of the qubit. However, 
it may be necessary to go beyond the dipole approxima- 
tion in the case of EWJN in a quantum dot. 



In this work, we present our study of the energy re- 
laxation of a single electron charge qubit in a DQD sys- 
tem and a single electron spin qubit in a single quantum 
dot. We assume that the primary source of field fluc- 
tuations are the metallic top gates of the quantum dot 
architecture. Back gates are typically a distance on the 
order of a micron from the qubits, which is too far to 
experience significant EWJN enhancement. We consider 
the detailed spatial variation of the electromagnetic field 
fluctuations and present results beyond the dipole ap- 
proximation which take into account the finite size of the 
quantum dot. We show that this extension of the dipole 
approximation removes the unphysical divergence in the 
relaxation rate at the metallic surface. This paper is or- 
ganized as follows: In Section |TT] we present our formal- 
ism for calculating the relaxation rate of a charge qubit. 
Results are presented for a DQD geometry. Section IIIII 
presents the formalism and results for the relaxation rate 
of a spin QD. In Section IIVI we derive an enhancement 
of the noise spectrum that results as the thickness of the 
metallic gate is decreased. Finally, Section |V] summa- 
rizes our results. Our results indicate that EWJN is the 
dominant cause of energy relaxation in some spin qubit 
experiments, particularly those performed in a small ex- 
ternal magnetic field, and is comparable in effect with 
previously studied noise sources in charge qubits. 
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II. CHARGE QUBIT 

We consider a charge qubit realized in a gated lateral 
DQD in an AlGaAs/GaAs hctcrostructure where elec- 
tron confinement along the z direction is much smaller 
than in the x ot y directions, so that we can safely de- 
couple the dynamics along x and y directions from the 
z direction^^. The total Hamiltonian of the charge qubit 
and its interaction with the electromagnetic environment 
is given by 



H — Hg + Hi, 



(1) 



where Hq is the Hamiltonian of the charge qubit in 
a DQD, which we model in the basis of the localized 
charge states {\L), \R)} as = e/2{\L){L\ - \R){R\) + 
A/2(|L)(i?| + \R){L\). e is the bias energy between the 
two dots, and A is the tunneling amplitude. In the en- 
ergy eigenbasis this Hamiltonian reduces to 



Ha 



(2) 



where az is the Pauli matrix, and hu = \J + A^. For all 
our calculations except those in Fig. [2l we will set e = 0. 
The interaction Hamiltonian i?i„t may be expressed in 
this same basis as 



following expression, which follows directly from Fermi's 
golden rule: 

(6) 

At finite temperature, the emission (transition from ex- 
cited to ground state) and absorption (transition from 
ground to excited state) rates are given by 

rl, = {l + N{oj,T))r,,, , Tl, =N{oj,T)T,,, (7) 

The Planck function iV(w,T) = l/[exp{hu}/kBT) - 1] 
gives the average occupation number of environment 
modes with frequency co at temperature T. The spec- 
tral density of the vector potential Xij is related to the 
retarded photon Green's function Dij byi^ 

X^, (r, r',oj) = JdT cxp-'"^([A, (r, t) , A, {r',t + r)]) 

= -— ImAj (r, f',w) , 
Co 

where i, j are Cartesian indices that run over x, y, z and 
the square brackets denote the commutator. Dij is ob- 
tained by solving 



H,nt Jdr [a. Air {r^ + <t. (rO] • i*(r , t) , (3) 

where A{r, t) is the vector potential of the fluctuating 
field. Mr and are associated with energy relaxation 
and pure dephasing in the charge qubit, respectively and 
are defined as 



Mr{r^=—i'l{f)ptlj^{r) 



ieh 



Mm 



mc 

e 



ip*f.{r) -0_ (r) V 



2mc 
ieh 



tpl (f) p ^+ (r) - tp*_ (f) p (f) 
ipl (f ) tp+ (r) V - i/'l (r) Tp- (r) V 



(4) 



Here m is the effective mass and p is the momentum op- 
erator of the trapped electron. Because we are operating 
within the weak field limit, the term proportional to 
has been dropped from the interaction Hamiltonian. We 
choose the gauge where the scalar potential = so 
that E = —dfA. The expression for Hint derives from an 
interaction in terms of operator quantities of the form 



-<5,, ( 



w^e (r, iu) 



did J 



{r,r' ,u) 



= -Airh 5^ {r-r')5^ 



jk- 



(8) 



Here the relative permittivity e(r, w) characterizes the 
geometry of a particular problem. In this section, we 
shall limit ourselves to the case where the metallic top 
gate of the lateral DQD is approximated by the half- 
space, z < 0. Then we can derive an analytical expression 
for D,J^ 



1 



27rift 

: f 



^\\Dij{k^ z, z' ^oj) dk , 



(9) 



Dx^{k,z,z' ,u) 



iq(z+z') 



7 7 

rs{k, uj) cos^ 9 —rp{k, ui) sin^ ( 



(10) 

where and rp are Fresnel's reflection coefficients given 
by 



Hint — 



2mc 



(l(r,i) ■p+p-A{r,t)'^ 



(5) 



This symmetrized version of the vector potential is not 
strictly necessary in our case since our qubit resides in the 
vacuum where \/ ■ A ~ 0, but we included it to keep our 
results more generally applicable. The zero temperature 
relaxation rate Fi^c = can be calculated using the 



rp{k,uj) 



eq-qi 



rs{k,u;) 



q-qi 



eq + qi q + qi 

Here k = (kx,ky), f|[ = (x — x',y — y'), 6 is the angle 
between k and the x-axis, and q = \J llP- j (p- — fc^ and 
qi — eijp- jf? — fc^ are the 2-components of the photon 
wavcvector in the vacuum and the metal, respectively. 
All other components of Dij can be derived from D^^)^ 
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In this work we consider e w icr /ixjeq, where cr = 6 x 10^ 
S/m is the conductivity of the copper gate. 

We pause briefly to mention that typical models of the 
interaction of a DQD with the electromagnetic field use 
the dipole interaction Hamiltonian 



H„,t = -E(r) ■ d, 
which will result in a relaxation rate of 



ri,c = 



(11) 
(12) 



in the quasistatic approximation, where d is the dipole 
moment of the qubit and E{r) is the strength of the 
fiuctuating electric field evaluated at the location of the 
qubit. This expression approximates that the electric 
field is uniform over the spatial extent of the qubit, which 
is equivalent to treating the qubit as a point dipole. As 
such, the qubit is able to couple to arbitrarily small wave- 
lengths of the electromagnetic spectrum, and the relax- 
ation rate is seen to diverge at shorter distances as ~ 1 / 
if the conductor is modeled with a local dielectric func- 
tioni^ii^. Using the complete electromagnetic interaction 
Hamiltonian ([S]) accounts for fluctuations of the field over 
the spatial extent of the qubit. If the wavelength of a 
particular Fourier component of the field fluctuations is 
smaller than the length of the qubit in that direction, its 
influence on the electron will average out and it will not 
contribute to qubit relaxation. The exact and dipole ap- 
proximation forms of the interaction Hamiltonian, Eqs. 
dS]) and (dD), converge when the distance from the gate 
becomes larger than the spatial extent of the qubit. 

We present calculations of the relaxation time for 
charge qubits that highlight the differences between these 
two forms of the interaction. First we consider DQDs in a 
one dimensional nanowire, which are realizable in semi- 
conducting nanowires^ or carbon nanotubes^. In such 
a geometry, the wave functions of trapped electrons in 
quantum dots have appreciable spatial extent in only one 
direction. We model the conflning potential of the DQD 
as a symmetric double square well potential and com- 
pute the lowest two eigenenergies and wavefunctions. We 
then compute the relaxation rate between these two low- 
est states which are separated by a fixed transition fre- 
quency w/27r = IGHz. A plot of the wave functions and 
the shape of the potential is shown in the inset of Fig. [1] 

We plot the energy relaxation time Ti vs. the distance z 
from the metallic gate to the DQD in Fig.[TJ In this plot, 
we choose the size of the dot in the x-direction d = 30 
nm and half the separation between the dots / = 30 nm 
(dash-dotted and dashed blue lines) and 60 nm (solid 
and dashed black lines). The curves that are shown in 
solid and dash-dotted lines are relaxation times for the 
exact form of the interaction Hamiltonian, whereas those 
shown in dashed lines are obtained using the dipole form 
of the interaction. The curves show significant deviation 
of the exact relaxation rate from the dipole relaxation 
rate at shorter distances and convergence of the two re- 
sults at longer distances. 
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FIG. 1: Energy relaxation time Ti vs. the distance from the 
metallic gate to the DQDs z for dot geometry of d = 30 nm 
and / — 30 nm (dashed and dash-dotted blue lines) and 60 
nm (solid and dashed black lines) at K, aj/27r = IGHz, 
and e = 0. Solid and dash-dotted lines are Ti times for the 
exact form of the interaction Hamiltonian, whereas dashed 
lines are for the dipole form of the interaction. The inset in 
the figure displays confining potential of a typical DQD in a 
one dimensional nanowire and the corresponding symmetric 
ground state and antisymmetric first excited state. 
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FIG. 2; Ratio of the Ti time for finite bias e to the Ti time 
at zero bias vs. e/A for three values of distances z from the 
metallic gate: 1 nm (solid blue line), 10 nm (dashed red line) 
and 50 nm (dash-dotted black line), for a dot geometry of 
d = 30 nm and I = 30 nm at K temperature. 



In Fig. [21 we present the ratio of Ti for a charge DQD 
qubit at bias e to the Ti obtained at e = versus the ratio 
e/A. An increase in bias increases the level splitting and 
decreases the dipole moment of the DQD. Since the relax- 
ation time Ti ~ l/ojcP, where d = 2/ sin (arctan (A/e)) 
is the dipole moment of the quantum dot, Ti increases 
for larger bias. 

Next, we present results from the relaxation rate cal- 
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FIG. 3: Energy relaxation time T\ vs. the distance from the 
metallic gate to the DQDs z for dot geometry of d = / = 30 
nm and I = 30 nm (dashed and dash-dotted blue lines) and 60 
nm (solid and dashed black lines) at K, ctj/27r = 1 GHz, and 
£ = 0. Solid and dash-dotted lines are T\ times for the exact 
form of the interaction Hamiltonian, whereas dashed lines are 
for the dipole form of the interaction. 



culation for a DQD in a two-dimensional quantum well. 
In this treatment we label the z-axis as the vertical con- 
finement direction and do not consider excitations along 
the z-direction. We model the confining potential by a 
symmetric double rectangular well in 2D and numeri- 
cally compute the lowest two eigenenergies and wave- 
functions. We then compute the electron relaxation rate 
between these two lowest states. The results are qualita- 
tively similar to the one-dimensional case and are shown 
in Fig. [21 where we plot the energy relaxation time T\ 
vs. the distance z from the metallic gate to the DQD. 
In this plot, we choose the size of the dot in the a;- 
direction d = 30 nm, the size in the y-direction / = 30 
nm and half the separation between the dots I to be 30 
nm (dashed and dash-dotted blue lines) and 60 nm (solid 
and dashed black lines). We find that for / = 30 nm 
and z = 90 nm, T\ is 4 [is while for I ~ 60 nm. Ti 
is 1.6 fis. These relaxation times are somewhat longer 
than the experimentally reported value of Ti = 20 ns 
in DQD-based charge qubits;^ We note that the relax- 
ation rate for a two-dimensional DQD is shorter than 
for a one-dimensional DQD of comparable geometry by 
about a factor of 5. A two-dimensional DQD is able to 
couple to obliquely oriented wavevectors in addition to 
those which point in the direction of separation between 
the dots, and this can be reasonably expected to enhance 
relaxation by a geometric factor of order unity. 



III. SPIN QUBIT 

We now focus on the calculation of the relaxation rate 
for a single electron in a spin qubit realized in a single 



quantum dot. Here the system Hg and the interaction 
Hint Hamiltonians are given by 



Hs = -gfiB S ■ Bo/2 
Hint = -g^J■B S ■ B{r, t) 



(13) 
(14) 



where a is the vector of Pauli matrices, g is the gyro- 
metric factor of the trapped electron in a quantum dot, 
/is is the Bohr magneton, Bq is the externally applied 
magnetic field and B{r, t) is the fluctuating EWJN field. 
The rate of spin flip from excited | f) to ground | 4-) at 
T = OK can be obtained from Fermi's golden rule 

Ti,. = ^ /"rfV / dV'A/,„,(f)A/,„,(r") 



X e,jkei.j'k'Xkk' {r, r',uj)njnj: 



(15) 



where repeated indices are summed over, and rij are the 
components of a unit vector h in the direction of Bq. The 
effect of finite temperature on the transition rates is the 
same as for charge qubits, as shown in Eq. ([7]). The 
magnetic spectral density xfj and the matrix element 



Mr,s{r) are 



SikmSpip Q^Q^ i:)„p(f, r", w) (16) 



(17) 



Here the spin qubit frequency uj = g^sBo/h and ipoi^ is 
the spatial part of the the ground state wave function of 
the spin qubit. Equation (|15p is a generalization beyond 
the dipole approximation of the simpler expressioni^: 



Ti, 



eoc^z 



(18) 



which has been obtained by using the quasistatic limit 
for the Green's function (jl6p . and assuming it is constant 
over the spatial extent of the qubit. Equation (fTS)) also 
assumes the external magnetic field Bq points in the z- 
direction. 

A plot of energy relaxation time Ti vs. the distance 
from the metallic gate z for a spin qubit is displayed in 
Fig. 21 Here we consider a single quantum dot of diam- 
eter c? = 60 nm and approximate the ground state spa- 
tial wave function of the spin qubit by the ground state 
wave function of a harmonic potential. We assume the 
Zeeman splitting between spin states is 50 GHz, typical 
of experiments in spin qubitsJi The solid line is the Ti 
time obtained using the non-local magnetic spectral den- 
sity while the dotted line is obtained for a local spectral 
density, which diverges as ^ 1/z as one approaches the 
metallic gate. The reason for saturation of the Ti time at 
smaller distances is similar to the case for charge qubits. 
There is a slight distinction in that the spin case involves 
a spatially extended dipole interaction, as opposed to 
the charge case which involves qenuine quadrupole and 
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FIG. 4: Energy relaxation time Ti vs. the distance from the 
metaUic gate to quantum dot z for spin qubit single dot ge- 
ometry of d = 60 nm at K temperature and g^BBo/2-nh = 
50 GHz. The solid line represent Ti time for the exact form 
of the interaction Hamiltonian, whereas the dashed line is for 
the dipole form of the interaction. 



far considered the simpler top gate geometry of a con- 
ducting half-space rather than the thin layer of finger 
gates used in these experiments. In the next section we 
address modifications to our calculations that we expect 
from more realistic gate geometry. 



IV. THIN METALLIC GATES 

A conducting half-space is an analytically convenient 
gate geometry, but a poor approximation to the thin top 
gates commonly used in semiconductor devices. In this 
section we present an exact treatment of the behavior of 
EWJN in the vicinity of a metallic film of finite thickness. 
Changing the half-space to a thin film affects EWJN by 
modifying the reflection coefficients and Vp. The power 
spectrum of the resultant EWJN is obtained by substitut- 
ing these modified reflection coefficients into the photon 
Green's function ([TU)) . and the relaxation time of, e.g. a 
charge qubit, is obtained by plugging Eq. pU)) into Eqs. 
© and The modified reflection coefficients for a film 
of thickness a take the form 



higher multipole contributions. This distinction is largely 
technical, however, and a saturation of Ti as z is ob- 
served in both cases. We find that the Ti time for a spin 
qubit in a GaAs quantum dot with an external magnetic 
field of 10 T and z = 90 nm is 150 ms which is larger than 
the experimentally reported valued of 0.55 ms, and gen- 
erally EWJN does not seem to be the dominant source of 
decoherence for semiconductor devices in large magnetic 
fields. GaAs has a strong spin-orbit interaction (SOI), 
which mixes the Zeeman-split spin states with orbitally 
excited states. Spin relaxation can then occur via cou- 
pling of the qubit to piezoelectric phonon noise in the 
2DEG layer. The relaxation rate from this mechanism 
scales as and is the dominant pathway for spin re- 
laxation at large external magnetic fields B > I7i2i^. 
Additionally, Marquardt and Abalmasso\^ calculate re- 
laxation of spin qubits from electric EWJN via the SOI. 
Again, mixing of the charge and spin states via the SOI 
allows spin relaxation to be induced from electric field 
fluctuations. They estimate the power spectrum of the 
Johnson noise using a lumped circuit model and found a 
B^ dependence of the relaxation rate. Our treatment in- 
volves a direct coupling of the fluctuating magnetic field 
from the top gates with the spin states, and our rate 
scales linearly with the magnetic field. We therefore ex- 
pect our relaxation pathway to dominate at low magnetic 
fields, and indeed while we predict a much slower relax- 
ation rate than measured by Amasha et ali^ for B ~ 7T, 
a.t B = IT our results predict Ti ^ 1.5 s which is compa- 
rable to their measured value of Ti = 1 s. Additionally, 
in Si quantum dots with a 2 T external magnetic field 
and z = 50 nm, we predict a Ti time of 6 ms which 
is smaller than the experimentally reported value of 40 
msJ^ However, it must be kept in mind that we have so 



gj + e-^q^ + 2iqqiecot(qia) 

11, ^ r - qj 

rs{k,u!,a) = 5 -. 

q + qf + 2iqqi cot(gia) 



(19) 



(20) 



They differ significantly from the half-space result only 
when the thickness a is of the order or smaller than the 
skin depth S, and they reduce to the half-space result for 
a ^ S. A derivation of Eqs. ((T^ and ([20]) is given in 
the Appendix. Equations and ([20]) are exact, but 
for a good conductor they can be cast into a simpler 
approximate form 



rp(k,uj,a) fa (iH — ^ cot (gic 
* ek 



rs{k,uj,a) 



1 - 



2c^qik 



cot((7ia) 



(21) 
(22) 



These expressions have been obtained by expanding 
Eqs. ([T9l) and (f20|) for large imaginary e and then tak- 
ing the quasistatic approximation q — > ik. The first ap- 
proximation is extremely accurate for copper near zero 
temperature and the second is accurate for all distances 
z such that EWJN is appreciably enhanced above black- 
body radiation^. The remarkable feature of Eqs. ((2T|) 
and (1221) is that they show the strength of the fluctu- 
ating fields outside the film are actually amplified rela- 
tive to the half-space result. This can be understood by 
analogy to the behavior of a particle trapped in a finite 
one-dimensional potential well. For a given width of the 
well, the wavefunction will have an exponentially decay- 
ing tail in the forbidden region. As the confinement is 
increased, the particle will be squeezed and its wavefunc- 
tion will leak farther into the forbidden region. It will be 
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FIG. 5: Ratio of energy relaxation time Ti from conducting 
film to Ti time from half-space vs. thickness of the film a 
for a DQD charge qubit in one dimension with dot geometry 
c? = 30 nm and / — 60 nm at K temperature. We take the 
exact form of the interaction Hamiltonian. The distance z 
from the film or half-space is chosen as follows: 2 = 10 nm 
(black dash-dotted line), z = 50 nm (blue dashed line) and 
z = 150 nm (solid red line). Other parameters are the same 
as in Fig. [1] 



interesting to see if this enhancement is observable in the 
Casimir attraction between 2 thin conducting plates. 

Using the modified expression for the reflection coeffi- 
cients, we compute the Ti time of a DQD charge qubit in 
one dimension due to the metallic film. In Fig. [51 we plot 
the ratio of the Ti time obtained for the film to the time 
computed for the metallic half-space as a function of the 
film thickness. We take the exact form of the interaction 
Hamiltonian for a variety of distances from the gate. We 
find that for distance z > a, the relaxation time due to 
the film can be reduced by over an order of magnitude 
relative to the half-space. It converges to the half-space 
result as z becomes smaller than the thickness of the film. 

Common semiconductor qubit architectures employ 
thin finger-shaped top gates which are more sparse than 
the films considered here. An exact treatment of EWJN 
from a detailed finger gate geometry would be pro- 
hibitively difficult, but we expect to a reasonable ap- 
proximation that EWJN from finger gates will be re- 
duced by a factor of the fraction of the top gate layer 
that is not composed of metal. Our results should then 
overestimate the relaxation rate by a geometric factor. 
We note however that newer accumulation-mode archi- 
tectures employ a second top gate above the confinement 
top gates2^. These accumulation gates are solid sheets 
and are typically around 100 nm from the qubit, so our 
treatment should accurately describe their contribution 
to relaxation. 



V. DISCUSSION AND CONCLUSIONS 

In conclusion, we have presented a detailed study of 
the effect of evanescent-wave Johnson noise on energy 
relaxation of quantum dots beyond the dipole approxi- 
mation. We have noted that previous studies of charge 
and spin qubits which use the dipole approximation al- 
low contribution from infinitely large components of the 
photon wavevector leading to overestimation and diver- 
gence of the energy relaxation rate as z — 0. We have 
demonstrated that it is possible to remedy this spurious 
divergence by taking into account the finite size of the 
quantum dot. While a non-local permittivity of the sur- 
face metal will remove the divergence in the field fluctua- 
tions at the boundary, we have shown that the finite size 
of the dot provides an alternative normalization mecha- 
nism by enforcing a finite cutoff in the magnitude of the 
contributing wavevector. In addition, we have derived a 
novel enhancement of the EWJN field fiuctuations that 
occurs outside a metallic film, relative to the field outside 
a metallic half-space. 

This manuscript has focused exclusively on relaxation, 
though we expect dephasing times from EWJN to be of 
comparable magnitude. The power spectrum of EWJN is 
linear in cj, and this will suppress contribution from the 
small frequency part of the electromagnetic spectrum, 
which typically enhances dephasing rates. While the 
temperature dependence of the relaxation rate is simply 
given by the Planck function, wc do expect a more non- 
trivial temperature dependence of the dephasing rate. 

Of particular interest are experimental signatures of 
EWJN-induced relaxation. Notably, at zero temperature 
the charge relaxation rate scales linearly with the qubit 
transition frequency and as the inverse cubic power of the 
distance between the qubit and the metallic top gates. 
The zero temperature spin relaxation rate scales linearly 
with the external magnetic field and inversely with the 
distance to the gates. 

Our results indicate that EWJN from the metallic top 
gate is not a dominant source of relaxation in charge 
qubits, but can be the dominant noise source for energy 
relaxation in spin qubits held at low external magnetic 
field. 



Acknowledgments 

We thank M.A. Eriksson for useful discussions. This 
work was supported by ARO and LPS grant no. 
W911NF-11-1-0030 and NSF grant DMR 0955500. 

APPENDIX 

Derivation of Green's tensor for a thin film 

Here we present the calculation for the retarded pho- 
ton Green's tensor outside of a thin conducting sheet of 
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permittivity e. The Green's function will satisfy Eq. ([8|). 
Here r' is simply a parameter for the purposes of solv- 
ing this set of equations, and wc take it to lie in the 
vacuum outside the conducting sheet. We will suppress 
the dependence of Dik{r,r',uj) on r' and uj to simplify 
the notation. The geometry of the problem is contained 
entirely in the permittivity function e(r, w). We take 
the boundaries of the conducting sheet to be located at 
z ~ —a and z = 0, with vacuum outside. Because the 
geometry is still translationally invariant in the x and y 
directions, we employ the same Fourier expansion Q as 
in Section |TT1 Solving Eq. ([5]) for a problem with planar 
symmetry is greatly simplified by separately considering 
the Fourier components of ([9]) that are polarized in the x 
and y directions. Dyy(f) may then be reconstructed as 



so that 



D 



(A-5) 

The term proportional to exjp{iq\z — z'\) is the free photon 
contribution to the power spectrum. It will have an imag- 
inary component and thus contribute to relaxation only 
in the radiative regime, k < co/c. Within a skin depth 
of separation from the metal, evanescent waves are or- 
ders of magnitude larger in field strength than these free 
photons. They may be safely ignored in this context. A 
similar calculation yields the result for Dyy,ky '■ 



Dyyir) = 



Ak-r\\ 



{^DyyM^ (k, Z) COS^ + Dyy,ky {k, z) Slu^ ( 



(A-1) 



Dyy,ky{k,z) = (r^(fc,c.,a)e*«(^+^') - e'^l^ 



where 



(A-6) 



where Dyy^k^ = Dyy when ky = 0, and Dyy^^y = 
Dyy when k^ = 0. The boundary value problem for 
Dyy^k^ (fc, z) then becomes 



Ae-'-?^ , z < -a 

iq 

(A-2) 



Our interest lies in the behavior of the fields for z > 0, 
so we need only to calculate C . Enforcing that Dyy^k^ and 
9Dyy^k^/dz are continuous across the boundaries results 
in 



C= ^r,(fc,w,a)e''9^ 



(A-3) 



where 
rs{k,uj,a) 



[q^ - ql) sin(gia) 



{ql + q^) sin(gfia) + 2iqqi cos(gia) 



2i singifl e 



q + n 



(A-4) 



rp(fc,w,a) 



e^q^ - ql) sin(9ia) 



ql + e^q^) sva{qia) + 2iqqie cos{qia) 



2i singifl e 



eq + qi 



eq-qi 



(A-7) 



A Taylor expansion of Eqs. (jA-4p and (|A-7P for large a 
in the evanescent range of wavevectors, i.e., a Taylor ex- 
pansion in powers of cxp(— 2|(7i|a), gives a monotonically 
increasing function of film thickness, a. However, a more 
careful treatment reveals that this is an error. The naive 
expansions of (jA-4p and (|A-7[) for large a neglect an en- 
hancement of the field spectrum that occurs for small k. 
In fact, EWJN is enhanced as the thickness is decreased 
for any good conductor. Specifically, the enhancement is 
preserved for a particular spatial Fourier component of 
the Green's function as long as |^| < 1. EWJN will 
eventually vanish as a — 0, but this does not occur until 
an unphysically small thickness is reached, on the order 
of 10"^'' m for copper at T = OK which is well below the 
applicability of the local permittivity model. 



^ D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

^ T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jeong, and 
Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 

^ J.R. Petta, A. C. Johnson, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Phys. Rev. Lett. 93, 186802 (2004). 

* J. Gorman, D. G. Hasko, and D. A. Williams, Phys. Rev. 
Lett. 95, 090502 (2005). 



^ Y. Hu, H. O. H. Churchill, D. J. Reilly, J. Xiang, C. M. 
Lieber, and C. M. Marcus, Nature Nanotechnology 2, 622 
(2007). 

N. Mason, M. J. Biercuk, and C. M. Marcus, Science 303, 
655 (2004). 

J.M. Taylor, J.R. Petta, A. C. Johnson, A. Yacoby, C. M. 
Marcus, and M. D. Lukin, Phys. Rev. B 76, 035315 (2007). 
* V. N. Stavrou and X. Hu, Phys. Rev. B 72, 075362 (2005). 



8 



^ M. Thorwart, J. Eckel, and E. R. Mucciolo, Phys. Rev. B 

72, 235320 (2005). 
° S. Vorojtsov, E. R. Mucciolo, and H. U. Baranger, Phys. 

Rev. B 71, 205322 (2005). 
^ D. C. B. Valente, E. R. Mucciolo, and F. K. Wilhelm, Phys. 

Rev. B 82, 125302 (2010). 
^ C. Henkel, S. Potting, and M. Wilkens, Applied Physics 

B: Lasers and Optics 69, 379 (1999). 
^ D. M. Harber, J. M. McGuirk, J. M. Obrecht, and E. A. 

Cornell, Journal of Low Temperature Physics 133, 229 

(2003). 

L. S. Langsjoen, A. Poudel, M. G. Vavilov, and R. Joynt, 
Phys. Rev. A 86, 010301(R) (2012). 
^ E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, 
Part 2, vol. 9 of Course in Theoretical Physics (Pergamon, 
1980). 

G. Agarwal, Phys. Rev. A 11, 253 (1975). 
^ J. M. Elzerman, R. Hanson, L. H. Willems van Beveren, 
B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwen- 



hoven. Nature 430, 431 (2004). 

M. Xiao, M. G. House, and H. W. Jiang, Phys. Rev. Lett. 
104, 096801 (2010). 

S. Amasha, K. MacLean, I. P. Radu, and D. M. Zumbuhl, 
and M. A. Kastner, and M. P. Hanson, and and A. C. Gos- 
sard, Phys. Rev. Lett. 100, 046803 (2008). 
'° A. V. Khaetsku and Y. V. Nazarov, Phys. Rev. B. 64, 
125316 (2001). 

F. Marquardt and V. A. Abalmassov, Phys. Rev. B. 71, 
165325 (2005). 

22 M. G. BorseUi, K. Eng, E. T. Croke, B. M. Maune, B. 
Huang, R. S. Ross, A. A. Kiselev, P. W. Deelman, L 
Alvarado-Rodriguez, A. E. Schmitz, M. Sokolich, K. S. 
Holabird, T. M. Hazard, M. F. Gyure, and A. T. Hunter 
Appl. Phys. Lett. 99, 063109 (2011). 

2^ In the ID case of a DQD realized in carbon nanotube^ or 
semiconducting nanowire,— we assume the electron dynam- 
ics along the x direction are decoupled from that of the y 
and z directions. 



